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ABSTRACT 
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00 : 

t:J" I We consider the constraints that can be derived from the spectral shape and 

^ I variability of TeV blazars on the homogeneous synchrotron self-Compton (SSC) 

<^ [ model. Assuming that the relativistic electron spectrum is a broken power 

i/^ [ law, where the break energy is a free parameter, we write analytical formulae 

^ I that allow to connect observable quantities to the physical parameters of the 

O ! model. We give approximate analytic formulae also for the case of Compton 

\ scattering occurring in the Klein Nishina regime, which is particularly relevant 

^ I for TeV blazars. In particular we find that even in the latter regime a power law 

Ph| component can be present at the highest energies. Further restrictions in the 

O I parameter space are set assuming that the break energy results from a balance 

\ between cooling and escape and that the soft photon lags measured in some 

I sources derive from radiative cooling of high energy particles. The constraints 

I can be summarized as allowed regions in the Doppler factor - magnetic field 

^ ' parameter space and are in principle sufficient to univocally determine the 

I model parameters and their uncertainties. 

We apply the method to three well studied sources Mkr 421, PKS 2155-304 
and Mkr 501. For Mrk 421 the available data are sufficient to fully constrain 
the model. The additional restrictions are found to be consistent with the 
model parameters, supporting the proposed interpretation. In the case of PKS 
2155-304, not yet detected in the TeV band, we estimate the peak frequency of 
the Compton component (~ 40 GeV) and the expected TeV flux. The derived 
physical parameters are similar for the two sources, with a relatively large value 
of the Doppler factor {S ~ 25) and a magnetic fleld of -B ~ 0.2 G. For Mrk 501 
we consider both the historical low state and the flaring state observed in April 
1997. In the flrst case consistency between the various assumptions is reached 
for 5 ~ 10 and B ~ 0.3 G. For the high state a similar value of the Doppler 
factor and a somewhat larger value of the magnetic fleld are indicated, while the 
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extremely large frequency of the observed synchrotron peak requires continuous 
injection/reacceleration of relativistic particles. 

Subject headings: BL Lacertae objects: general-gamma rays: observations- 
gamma rays: theory- radiation mechanisms: non-thermal-X-rays: theory 

1. Introduction 

The discovery of intense gamma-ray emission from blazars had a strong impact on 
our understanding of relativistic jets. Independent evidence of relativistic beaming was 
provided by the condition that the source should be transparent to high energy gamma-rays 
(Mattox et al. 1993, von Montigny et al. 1995, Dondi & GhiseUini 1995). The gamma-ray 
fluxes and spectra, together with data at lower frequency allowed detailed estimates of the 
physical parameters of different models (Maraschi, GhiseUini and Celotti 1992; Blandford 
1993; Sikora, Begelman and Rees 1994; Dermer and Schlickeiser 1993). 

More recently the discovery of TeV emission from some objects (Mrk 421, Mrk 501 and 
lES 2344-1-514) and the detection of very fast variability at these high energies motivated 
new theoretical and observational studies (e.g. Celotti, Fabian and Rees 1998) 

The objects detected at TeV energies have similar spectral energy distributions (SED) 
from radio to X-rays. They emit synchrotron radiation with power peaking in the X-ray 
band. For Mrk 421 and Mrk 501 the X-ray and TeV radiation vary in a correlated fashion 
(Macomb et al. 1995, Buckley et al. 1996, Catanese et al. 1997, Schubnell 1997) suggesting 
that the TeV emission derives from the same high energy electrons responsible for the 
synchrotron X-rays. For Mrk 421 the TeV flux has been observed to vary with a 20 
min timescale (Gaidos et al. 1996), which strongly limits the size of the emitting region. 
Furthermore a lag has been measured between soft and medium energy X-ray photons. A 
lag of the same order between soft and medium X-ray photons has been measured also in 
PKS 2155-304, whose spectral energy distribution closely resembles that of Mrk 421 up 
to the X-ray range, but has not yet been observed at TeV energies due to its southern 
declination. The lags strongly suggest a radiative origin yielding a new important constraint 
on the physical parameters of the emitting region. For Mrk 501, similar in spectral shape 
to the previous two, an extraordinary flaring behaviour has been observed recently (Plan et 
al. 1998). 

The simplest model proposed for TeV blazars involves a single zone and a single 
population of relativistic electrons emitting synchrotron radiation from radio to X-rays 
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and inverse Compton radiation from X-rays to gamma-rays. For this class of sources the 
seed photons for the inverse Compton process are hkely to be the synchrotron photons 
themselves (SSC) (see e.g. GhiseUini & Maraschi 1996, Mastichiadis & Kirk 1997), although 
other models have also been considered (for a review see Sikora 1997). We will limit 
ourselves to the SSC model in the following. 

Our scope is to provide general analytic expressions for all the constraints that can be 
envisaged within this simple model and discuss the consequent restrictions in the parameter 
space for the above mentioned sources individually and as a group. A similar approach 
was initiated by Ghisellini, Maraschi and Dondi (1996) (hereafter GMD), however it did 
not include all the conditions and was limited to regimes well below the Klein Nishina 
limit, while for TeV blazars the Klein Nishina corrections are important. Work along these 
hues has been developed by Dermer, Sturner and Schlickeiser (1997) and by Bednarek 
and Protheroe (1997). The latter authors chose to use precise formulas requiring however 
numerical computations, with loss of generality. 

The paper is organized as follows. In section 2 we discuss the constraints that can be 
put on the homogeneous SSC model using all the available observations. In section 3 the 
relevant formulae are derived which link model parameters to observative quantities. The 
effects of the Klein Nishina limit are discussed in detail and revised analytic formulae are 
given valid also in the latter regime. In section 4 we apply our diagnostics to the case of 
three well known blazars and finally (Section 5) we discuss the results. 

2. Model constraints 

The homogeneous SSC model assumes that radiation is produced in a single zone of 
the jet (approximated here for simplicity as a sphere with radius R), relativistically moving 
at small angle 6 to the observer's line of sight. Photons up to the X-ray range are produced 
by relativistic electrons through the synchrotron process and are subsequently inverse 
Compton scattered by the same electrons to energies in the 7-ray range. The observed 
radiation will be strongly affected by relativistic effects. The key parameter is the Doppler 
factor 5 = [r(l — f3 cos^^)]~^, where P is the bulk Lorentz factor and (3 = v/c. 

The observed SED of TeV emitting blazars (in a vF^^ plot) is characterized by two 
broad peaks, the first approximately in the X-ray band, the second in the 100 GeV-1 TeV 
band . Below and above the peaks the spectrum is quite smooth and can be approximated 
with power-law profiles with indices q;i,q;2 respectively where cci < 1 and 0:2 > 1 (see e.g. 
GMD). Note that within the present uncertainties the "left" side of both peaks can be 
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described by the same spectral index. Note also that the spectrum in the TeV range is 
poorly known and may be affected by intergalactic absorption (see e.g. Stecker and De 
Jager 1997 and reference therein). In figure (|I|) we show the observed SED of Mrk 421 at 
two different epochs (Macomb et al. 1995). 

The observed spectral shape requires that the the relativistic electron spectrum 
steepens with increasing energy. We approximate this behaviour with a broken power 
law (GMD) with indices ni < 3,n2 > 3, respectively below and above the "break" energy 



With these approximations we can completely specify the model using 7 parameters: the 
magnetic field intensity B, R, 6, the slopes ni and n2, the Lorentz factor of the electrons 
at the "break", 7;,, and the electron density parameter K. The peak synchrotron power 
is emitted by electrons with the break energy Ef, = 7bmc^. The maximum energy 7maa: 
attained by electrons as well as a possible lower limit 7mm are unimportant here provided 

that 7max > lb and Jrnin < 100 

Spectra computed with the above assumptions are compared with the Mkn 421 data 
in Fig. 1. The theoretical curves show that the double power-law approximation can 
reproduce reasonably well the shape of the observed SED. 

The available data on the SED can be used to derive 6 "observable" quantities of 
particular relevance for the above model, namely the indices ai and 02, the frequency of the 
synchrotron and inverse Compton peaks, and u^, and the peak luminosities Ls{vs) and 
Lc(z/c). These quantities are rather uncertain being ideally the result of a fitting procedure 
on simultaneous multifrequency data, but in practice often inferred from observations 
obtained at different epochs and with incomplete frequency coverage. The best sampled 
spectral region is usually the optical to X-ray range, therefore ai, a2, I's and Ls{h's) are 
better determined than z/^ and Lcii^c)- Fig. 1 can serve as an illustration of the uncertainties 
in the determination of the " observables" for Mrk 421, perhaps the object with the best 
overall data. We recall that the one-zone homogeneous model is selfabsorbed at radio 
frequencies and cannot explain the radio emission which implies further contributions from 
the outer regions of the jet. 

It is important that, at least in principle, one can derive from the broad band spectral 
observations 6 quantities that have to be reproduced with the 7 parameter model. One more 
observable quantity is sufficient to close the system. This can be provided by the minimum 
timescale of variation, t^ar, which can be directly connected to the source dimension, R, 
through the causality relation R < ct^ar^- This relation gives only an upper limit on 
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the dimension of the source, but a small dimension implies a high Doppler factor if the 
transparency condition is taken into account (see below). In order not to force S to be too 
large we will adopt limits to the radius of the emission region in the range ct,ua,.5- ^ctyar^- 
At this point the system is closed and any additional constraint that can be satisfied should 
be regarded as supporting the model. These are: 

• the transparency of the source to 7-rays. 

The high energy photons may interact with the low energy ones producing pairs. The 
fact that we observe the TeV emission poses a strong limit to the optical depth of the 
source to photon photon interaction and therefore to the energy density of the soft 
radiation. Prom this limit one can obtain a lower limit on the value of the Doppler 
factor (see e.g. Dondi & GhiseUini 1995). 

• the radiative interpretation of the time lags in the light curves at different frequencies. 

Recently, observations of Mrk 421 and PKS 2155-304 with ASCA (0.3-10 keV) have 
revealed that when pronounced flares occur, the soft photons lag the harder ones ( 
Takahashi et al. 1996, Urry et al. 1997). The simplest way to explain the time lags in 
a homogeneous model is through the hypothesis that they are associated to the time 
necessary for freshly injected high energy electrons to cool. With this assumption the 
measured time lags determine the cooling time setting a well defined condition on the 
model. 

• the consistency of the break energy of the electrons with the condition of balance 
between cooling and escape. 

In our discussion the electron energy distribution is not calculated self-consistently 
(for an example of this alternative approach see e.g. Inoue and Takahara 1996, 
Mastichiadis and Kirk 1997, GhiseUini et al. 1998). The two spectral indices and 
the break energy arc determined from the observed spectra. Having determined at 
the same time the dimension of the emission region and the magnetic and radiative 
energy densities we can verify a posteriori whether the cooling time of the electrons 
at the break is of the order of the escape (expansion) time from the emission region. 
We will assume that consistency is achieved if the required escape/expansion velocity 
is in the range c/3 — c. We recall however that in the standard hypothesis of a power 
law injection and subsequent cooling and escape in a homogeneous region, a spectral 
change Aa — 0.5 is expected, while in general for the sources discussed below the 
observed steepening is larger, suggesting a more complex situation. 

We approximate the cooling time as the minimum between the synchrotron and the 
Compton cooling times. In the B-S parameter space the regions of validity of these 
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two approximations are separated by the line expressing equipartition between the 
magnetic and radiation energy density. In each of these regions a value of the cooling 
time defines only one line. Bednarek & Protheroe (1997) discuss a similar constraint 
on the cooling time of the electrons, but they incorrectly apply the condition to 
both the synchrotron and the Compton cooling times simultaneously obtaining two 
independent constraints. 

Summarizing, the SSC homogeneous model needs 7 independent parameters: the spectral 
observations give 6 independent quantities, and the seventh observational information used 
to close the system is the minimum variability timescale. Other constraints are provided by 
the optical depth for 7-7 interaction, by the observed time lags between soft and hard X-ray 
photons and by the consistency of the value of 75. These last conditions are additional 
independent constraints: if verified, these constraints argue for the physical consistency of 
a homogeneous SSC model. 



3. Relations between the model parameters and observables 

In this section we derive the fundamental relations between the model parameters and 
the observable quantities discussed above. We chose to do this analytically at the cost of 
some inevitable approximation in order to maintain generality. Most of these formulae are 
also discussed in GMD and Dermer ct al. (1997). 

We assume that the two spectral indices can be determined directly from the data, 
independently of other quantities. The structure of the relations is then such that we can 
eliminate 3 of the remaining 5 parameters (i.e. 7;,, K and R) and express the constraints as 
allowed regions in the B-6 parameter space. 

3.1. Peak frequencies 

In the scheme that we consider the dominant synchrotron power is emitted by electrons 
at the break (those with Lorentz factor 7^). Primed quantities refer to the blob frame so 
that the observed frequency is given by i>s = i^'s^- The photons at i>s dominate the energy 
density of the target photons and the inverse Compton power in the Thomson limit will be 
dominated by photons at frequency Uc — (4/3)7^i/^5, produced in the scattering between 
the electrons at the break and the photons with frequency u'^ (see also GMD). Inverting the 
last equation one obtains: 
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Therefore the ratio between the peak frequencies gives directly the value of 7;,. 

The synchrotron frequency averaged over the spectral shape for an electron of Lorentz 
factor 7fe is: 

i/, = 3.7- 10^25^ (3) 



using the value of 7^ given by eq. (0), gives: 



B6 = (l + z) (4) 

The last formula shows that for fixed and u^, B and 5 are inversely proportional. We will 
see in the next section that in the KN limit this situation is reversed. 



3.2. Peak luminosities 

Another relation can be obtained from the ratio of the total luminosity of the 
synchrotron peak and the total luminosity of the self-Compton peak, which is directly 
related to the ratio between the radiation and the magnetic field energy densities inside the 
source (in the comoving frame): 

T W 

For our aims it is useful to express the total luminosity through the peak luminosity, which 
is more accessible to observations. A simple calculation shows that the total luminosities 
{Lg and L^) and the peak luminosities [Ls{h's) and Lc^Uc) ] are related by: 



Ls,c = /(«!, <y2)^s,cLs,c{'^s,c) (6) 



where /(ai,Q;2) is given by 
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f{ai,a2) = -^ — + — ^ (7) 
1 — «! a2 — i- 



Expressing U'gy^ as: 



UL 



we can rewrite eq. (Bl) as: 



If we express the source dimension R with the variabihty timescale 



R < Ctt,ar(5(l + z) 



(10) 



we can finally write: 



>{l + z) 



2(z/,L,(t/,))V(ai,a2) 



1/2 



This formula relates the magnetic field intensity and the Doppler factor to observable 
quantities. Eq. ([T^) is an inequality because of the use of the inequality eq. (PUj); taking 
ctvar/^ as a lower limit on the radius of the source, as discussed in Section 2, eq. (|rT|) will 
be represented in the -B, 5 space by a strip of finite width. Thus the intersection of eq. flTlJ 
with eq. allows to estimate B and 5. 



4. The Klein Nishina regime 

4.1. Spectral shape of the Inverse Compton emission 

The adopted spectral shape for the electrons ensures that the synchrotron spectrum 
has a well defined peak (break) (as observed). Thus both, electrons and photons, have a low 
energy branch and a high energy branch (below and above the break). With the adopted 
approximations one can compute the inverse Compton spectrum as the sum of 4 separate 
integrals: /i,i,/i,2 and J2,i, /2,2, representing the contributions from electrons of the low 



- 9 - 



(1) and high (2) energy branches scattering photons of the low and high energy branches 
respectively. Full expressions are given in the Appendix. This subdivision is very useful in 
the discussion of the Klein Nishina limit. 

The different contributions, computed using the well known 5-function approximation 
for the single electron emission and a step function for the energy dependence of the cross 
section, a = cr^ for 7Z/' < 3/4 and cr = otherwise (see Coppi and Blandford 1990 for 
a critical discussion), are shown in Fig. 2. All the spectra are computed in the source 
frame for fixed values of u'^ = 10^^ Hz and varying 7?, from 10^ to 10®, that is for increasing 
importance of the KN limit. 

Since the assumed particle spectra have no high energy cut offs, at high enough energy, 
the KN limit is always relevant and steepens the spectrum. The limiting frequency can 
be estimated considering that the largest contribution to the radiation energy density 
is due to photons of frequency I'g. Therefore the KN limit will set in at a frequency 
z/^ > {mc^/h)^/i>'g. This is seen in Fig. 2a and 2b, where three (12, 21, 22) of the four 
components in the total spectrum cut off at high energies, while the peak of the inverse 
Compton spectrum occurs as expected in the Thomson limit. 

However the effect will be severe and affect both the value of the peak frequency and 
the bulk of the emitted power if uk is close to uq, that is if the KN regime sets in already 
for the scattering of electrons of energy 7^ with photons at the peak of the synchrotron 
emission. Therefore the condition for the KN regime to affect severely the luminosity and 
the peak frequency of the inverse Compton emission, requiring a revision of the formulae 
given in section 3 is: 



, 3 mc^ , , 

> J— (12) 

Using eq. (@) for to express the above condition with observed quantities, eq. (|l^) 
becomes a condition on 6: 



6 <6 



KN 



(3/4)(mc2//i)^ 



1/2 



(13) 



Since eq. (0) (valid in the Thomson limit) gives the lowest possible value of 7?,, eq. (|13|) 
gives the limit above which the full KN regime sets in. However even for 5 somewhat above 
the KN limit corrections are non negligible (see Section 4.2) 



The contribution involving both electrons and photons on the high energy branches 
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( /2,2 ) is already strongly suppressed in Fig. 2b and disappears completely when the 
KN condition is met (Fig. 2c, 2d). The other contributions are reduced according to the 
availability of photons below the KN threshold. In Fig. 2c and 2d /i^2 is negligible and the 
low energy side is dominated by by the Ji.i branch. 

It is interesting to note that for the assumed spectral shape the scattering of highest 
energy electrons with very low energy photons gives a significant power law contribution 
to the inverse Compton spectrum at the highest energies. The latter component is not 
important energetically (at the peak) however it provides a high energy tail which extends 
the range of the inverse Compton emission. The spectral index of the latter contribution can 
be easily derived from the expression for J2,i given in the Appendix and is axN = 2^2 — «i- 
This component could explain the power law spectra measured in the TeV range for Mrk 
421 and Mrk 501 (Zweerink et al. 1997). 



4.2. Compton peak frequency and luminosity in the KN limit 

In the KN regime both the expression of the Compton peak frequency and luminosity 
must be modified. As shown in Fig. 2c, 2d, with the assumed electron spectrum the peak 
frequency in the KN limit is determined by the scattering between all the electrons and 
photons with v < Vg. Maximizing the sum of Ji^i and l2,i (see Appendix) leads to the 
following relation between and 7^: 



, 5 mc^ , 5 



'1 + z h 



l + z 



where: 

5f(ai,a2) = exp 
is a factor smaller than one. 



1 1 

+ 



«! — 1 2(q!2 — «l) 

Substituting 7;, derived from eq. (IT^ in eq. (H) we obtain: 



(14) 
(15) 



5 \ h 3.7 -106 l + z 



(16) 



which is the analogous of the eq. (^) for the KN regime. Note that in this equation B and 
S are directly proportional. 



Having obtained a revised expression of eq. (|^) we can use the latter (eq. [|T^) to 
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express (|T^) through observables. We therefore arrive at a revised version of the threshold 
for the KN regime: 



6 < 5kn 



(3/4)(mc7/i)2 



1/2 



'1/2 



(17) 



Since 5'(ai,a2) is less than 1, eq. ([T7| ) is more stringent than eq. (|13]): therefore we 
find two different limits dividing the regions of applicability of the two treatments. Between 
this two limits there is a transition region, where both the approaches are approximations. 

In the KN regime the power in the Compton peak is substantially diminished, equation 
(H) should be replaced by: 



TJl 

syn,avail 



(18) 



where UL 



syn, avail 



includes only photons below the KN limit for electrons with the break 



energy, that is with v'^ < mc //17b. We can write: 



syn, avail 



3r?ic^/4/i7i, 



syn 



(19) 



where e^y„(z/Q) is the energy density of the synchrotron radiation for a given frequency. A 
simple calculation shows that the integral is given by: 



n' 

^ syn, avail 



(20) 



where the factor in parentheses represents the ratio between the available energy density 
and the total energy density. From the last relation we can obtain the following expression 
for the energy density ratio in the comoving frame in the KN regime: 



U' 

^ syn 



hi 



l^sK: 



g[oii,a2) 
l + z 



oi — 1 



(21) 



Using the last formula and expressing the total luminosities through eq. (H) and the radius 
with eq. (|10[) we can finally obtain the analogous of eq. (0) for the KN limit: 
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(l-ai)/2 - 


g{ai,a2) 







2(^/,L,(z/,))V(ai,a2) 



1/2 



3mc" 



1— ai 



(22) 



5. Additional constraints on the model parameters 

5.1. The pair production opacity 

Another strong constraint is obtained from the condition of transparency of 7-rays to 
pair-production absorption. Following Dondi & Ghisellini (1995) we can write this condition 
in the form: 



s > 



ax 



dl{l + z) 



2/3 



1/(4+2/3) 



(23) 



where z/q = 1.6 • 10"^^ /u^ is the frequency of target photons and /? is the spectral index of 
the target photons {ai for i/q < and a2 for > Us). The transparency condition does 
not depend on the specific emission mechanisms: therefore it provides a strong independent 
constraint on the minimum value of the Doppler factor in the SSC model. 



5.2. Cooling time 

The radiative cooling time for a synchrotron-self Compton emitting electron is given 



by: 



cool 



4 arc 
3 mc^ 



1{U'b + U'syn,avail) 



(24) 



Clearly the cooling time is determined by the fastest of the two cooling processes that is 



synchrotron cooling for U'^ > Ui.. 



syn,avail 



or inverse Compton cooling for U'^ < U'g. 



syn,avail • 



Let 



us first consider for simplicity the Thomson regime. In the case in which the magnetic field 
energy density dominate we can approximate the cooling time as (after some manipulations): 



(25) 
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cool 



comp 



1/2^13/2^-1/2 



(26) 



where Cs and Cc are constants. 



In the observer frame: 



t 



sync 



^^^-3/2^-1/2^-1/2. 



comp 



(27) 



The equipartition condition U'^ = U'^^^ divides the regions of vahdity of the two 
approximations in the B — 5 plane. The Compton coohng approximation will be valid in the 
lower region (where U'^ < U'^y^^), while the synchrotron cooling approximation will be valid 
in the upper region. In both cases a given value of tcooi defines an inverse relation between 
B and 6, that is lines with negative slopes, intersecting each other at the equipartition line 
(see Figures in the Application sections). Thus a condition on the cooling time corresponds 
to a minimum value of 6 above which the condition is represented by a double valued 
function with two branches. Lower values of 6 correspond to shorter cooling times. 



The measurement of time delays between variations in different energy bands provides 
an important restriction on the possible values of the physical parameters of the model. 
The most natural way to explain the lags in a homogeneous model is to interpret them 
as due to the cooling time of the emitting electrons (Takahashi et al. 1996, Urry et al. 
1997). In this picture the quiescent emission is produced by the "stationary" population 
of electrons, while the flare is due to the injection of a monoenergetic population of high 
energy electrons: while these electrons cool their emission drifts to lower frequencies as 
observed. Adopting this hypothesis we can relate the observed time lags to the cooling 
time of the electrons and therefore to the physical parameter of the source. If the cooling is 
dominated by the synchrotron emission (i.e. U'^ > U'^y^, where U'^y^ is U'^yn^avaii 
limit), we obtain the following relation (see also Takahashi et al. 1996): 



5.3. Time lags 




1/3 



2/3 



G 
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where and z/i are the X-ray frequencies { uq > ui) and Tobs is the measured lag (in 
seconds). In the Compton coohng dominated region ( i.e. U'^ < U'^y^ we can write the 
equivalent relation: 



= ^^i^ 

l + z 



Tobs 



1 2 



„3+2 
^ ''var 



1 - (^1/^0)^/2 



(29) 



5.4. Cooling vs. escape 

In our approach the value of 7^ enters as a free parameter in the model and is 
determined from the observations. However in more sophisticated models 7b could be 
determined selfconsistently by the equilibrium between injection, cooling and escape of 
electrons from the source. Here we include a posteriori the condition that the cooling time 
at 7f, is equal to the escape time. The latter constraint can be written as: 



R 



4 c^tc , , 
3 mc^ 



U' \ 

^ syn, avail) 



5 • 10^ 



syn,avail 



(30) 



B) 



where Peso is the electron escape velocity in unit of c and U'^ and Ui.yn^avaii ^^6, respectively, 
the magnetic field energy density and the radiation energy density. 
If cooling is dominated by the synchrotron process, expressing (|30|) in terms of the 
observational quantities and Uc gives: 



> 



1/2 



1/2 

esc 



(31) 



for the Thomson limit and: 



B > 



5 • 10^5'(q;i, 0^2) 'mc^ 1 



for the KN limit. 



1/2 



A 



1/2 

esc 



(32) 



Even when Compton cooling is less important than synchrotron cooling, if it occurs in the 
KN regime, the value of 7f, must be taken from eq. (|1^) instead of eq. (|^), which causes the 
difference between (|3Tf) and (B^). 



In the Compton-cooling dominated region the equations are: 
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H 2.1 



B6'^ > 5.2 ■ 10 



-14 




8/5 



-13/5 



(33) 



var 



for the Thomson regime and: 



-|l/(l-ai) 



=e(«l,«2) 



K(l + z)] 



(3ai-2)/(l-Qi) 



(34) 



[^(ai, 02) is a constant given in the Appendix] for the KN regime. 

The electron escape time is not a well known parameter and depends on the transport 
processes in the emitting region: in the applications we check the consistency of the 
observed value of 7{, for values of 13 esc in the range 1 — 1/3. 



In the previous section we have shown that the constraints considered lead to equations 
that can be expressed in terms of two parameters: the magnetic field and the Doppler 
factor. Each constraint then corresponds to an allowed region or to a line in the Log B - 
Log 5 plane. Three equations do not apper explicitly: one determines 7^ from the peak 
frequencies, the second determines the radius from the variability timescale, the third the 
normalization constant of the electron spectrum from the observed synchrotron luminosity. 
We examine here what consequences can be derived comparing the constraints available for 
the three brightest and best observed BL Lacs with synchrotron peaks near or within the 
X-ray band. 



We start from the following values for the observable quantities, derived from the 
simultaneous observations of Macomb et al. (1995) (the corresponding data are shown in 
Fig. 1) : z/, = 3 • 10^^ Hz, = 6.5 • 10^^ Hz, = 0.5, = 1.75 i/sL,(z/,) = 4.7 • 10^^ erg 
s~^ and VcLc{v^ = 6.5 • lO^'' erg s~^ (see Table 1). These refer to the low state (for the 
high state we find similar results for the values of B and 6). The TeV observations show 
extremely short variability timescales (Gaidos et al. 1996 report a flare with tyar — 20 min), 
while in the X-ray band variability timescales of several hours are observed (see Takahashi 
et al. 1996): in Fig. (|^) we have used t^ar — 1 h. For Mrk 421 the soft photon lag within 



6. Applications 
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the ASCA band was measured by Takahashi et al (1996) yielding Tabs = 3200 s between 
Ei = l keV and E2 = 5 keV. 

In Fig. @ we show the parameter space for Mrk 421. The transparency condition 
(^) calculated for photons with u = 10^^ Hz yields a lower limit 6 > 15 represented by 



the vertical line. The KN condition [eq. ([T3|) ] yields an upper limit 6 < 14, while the 
more stringent version [eq. ([l^)] yields 6 < 31. For 14 < 5 < 31 we are therefore in an 
intermediate region where both limits are inaccurate. We chose to use the Thompson limit 
here since the resulting value of S will be 25, close to the upper end. We recall however that 
even if the Thomson limit is adequate for the peak of the Compton component, KN effects 
cannot be neglected in the computation of spectra at higher energies (see Fig. 1, Fig. 2a). 

The observed ratio between the two peak frequencies [eq. (^] imposes an inverse 
proportionality relation between B and 6. This is shown as a shaded area allowing for a 
factor 3 uncertainty in the actual value of the Compton peak frequency. 

The ratio of the peak luminosities (|TT]) yields a steeper relation between B and 6, 
B oc 5^^. Again it is represented by a shaded area to allow for uncertainties in the 
observational input data. 

The intersection between the two regions above defines the allowed region, where the 
actual values of the parameters should fall. This region could shrink if the observational 
data improve, but given the present gap between the energy ranges of EGRET and the 
ground based Cherenkov telescopes it is unlikely that it could be substantially reduced. 

The time lag constraint has a shallow dependence on 5 ( i? oc 5^^^^) and is well defined 
observationally. Therefore it is represented as single thick line (labelled / ) in Fig. 3. It 
is extremely interesting that this line, which is completely independent of the spectral 
constraints, is consistent with the parameter region allowed by those constraints. It is also 
interesting that the region defined by the condition that the break energy of the electrons is 
determined by a balance between cooling and escape eq. ( |3TD (region C) encompasses the 
lag constraint in the relevant parameter interval. 

We conclude that for Mrk 421 the spectral constraints are consistent both with the 
radiative interpretation of the soft photon lag and with the break energy as due to cooling. 
These are persuasive arguments in support of the simple general model. The estimated 
parameter values are B ~ 0.25 and 5 ~ 25 (see Table 2) 

Note that the value of 5 is rather high compared to other estimates (Ghisellini et al. 
1998) but similar to that found from direct model fits by Mastichiadis & Kirk (1997) and 
Ghisellini & Maraschi (1996). 
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6.2. PKS 2155-304 

This BL Lac object has been studied for a long time with numerous multifrequency 
campaigns (Urry et ah 1997; Edelson et al. 1995). It has been detected as a 7-ray source by 
EGRET on one occasion (Vestrand et ah 1995) and more recently in a flaring state (lAU 
Circ. 6774). No simultaneous X-ray and gamma-ray observations are available yet. It has 
not been detected in the TeV band. This may however be due to its southern declination, 
since in the southern hemisphere Cherenkov telescopes have started observations only 
recently. 

The parameters of the synchrotron component are inferred from the spectrum measured 
during the 1994 multifrequency campaign (z/^ = 10^^ Hz and VsL{vs) = 1.5 ■ 10^^ erg s~^, 
Urry et al. 1997). Comparison with previous observations (see Sambruna et al. 1994) shows 
that during these observations PKS 2155-304 was in an intermediate state, a factor 2 less 
luminous than the brightest state recorded by EXOSAT. 

The gamma-ray spectrum measured with EGRET is hard, implying that the Compton 
peak is beyond the EGRET range. However the absence of TeV observations does not 
allow to determine the position of the peak. We will therefore use the other constraints to 
estimate the peak frequency of the Compton emission and the associated TeV luminosity, 
assuming that all of them should be satisfied as is the case for Mkn 421. 

In Fig. (^) we show the parameter obtained for Uc = 10^^ Hz and for an "observed" 
luminosity of VcLivc) = 1-5 • 10^^ erg s~^ extrapolated from the GeV range. We assume 
a minumum variability time of t^ar = 2 • 10^ s and a delay of Tobs — 5 • lO'^s between the 
variations observed in the 0.75 keV and 6 keV energy bands (Urry et al. 1997). The vertical 
line representing the transparency condition is plotted for photons with frequency v ~ 10^^ 
Hz for consistency with the following estimates of the TeV flux. Note, however, that the 
lower limit on the Doppler factor given by the transparency condition would be weaker in 
the absence of TeV emission. 

In this case, the largest uncertainty is on the position of the Compton peak, while 
the frequency of the synchrotron peak is reasonably well known. We therefore plot the 
uncertainty band for fixed and for ranging from 2.5 • 10^^ Hz to 2.5 • 10^^ Hz. For these 
values of the peak frequencies the KN limit yields 5kn < 3, smaller than the transparency 
limit: we can make use of the Thomson formulation. 

We find that all the constraints are satisfied for values of B and 6 very similar to those 
found in Mrk 421 {B ~ 0.25 and 6 ~ 25); the cooling-escape consistency for 7^ is satisfied 
also in this case. 
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We can estimate the value of the TeV flux of PKS 2155-304 with the hypothesis that 
beyond the break the spectral indices of Mrk 421 and PKS 2155-304 are similar. From the 
spectrum of Mrk 421 we can obtain the value arev — 2.2; from the inferred value of the 7 
peak flux of PKS 2155-304 we can evaluate the TeV flux as: 



that gives an integrated flux F{> ITeV) ~ 10~" ph s~^ cm ^, about 25% that of Mrk 
421 in the in the non-flaring state. Hence we expect that a detection of PKS 2155-304 
should be possible but requires sufficient sensitivity, due to the low flux. TeV observations 
of PKS 2155-304 would be extremely important also because they may reveal or constrain 
absorption of high energy photons by the intergalactic IR background (Stecker and De 
Jager 1997), better than in the case of Mrk 421 because of the larger distance. 



The BL Lac object Mrk 501 (z=0.034) underwent a dramatic TeV and X-ray outburst 
in April 1997 (see e.g. Plan et al. 1998, Catanese et al. 1997). Previous observations 
showed a spectral distribution very similar to Mrk 421, with a synchrotron peak in the soft 
X-ray region; in contrast, during the April observation the peak was found at about 100 
keV, a shift in energy of at least a factor 100. The broad band spectral shape leads to 
interpret the hard X-ray emission as synchrotron radiation and the flare as due to injection 
of electrons with very high energy (Plan et al. 1998). 

In Fig. (^) we show the parameter space for Mrk 501 in the quiescent state. In this 
state we can apply the relation for the Thomson case, given that the value of S required by 
the transparency condition is higher than 6k n- We have used the following values for the 
observable quantities: = 10^^ Hz, z/^ = 10^^ Hz, VgLs^Us) = 10'^'^ erg/s, UcLdi^c) = 5 • 10^^ 
erg/s, tyar — 10 h, ai = 0.5 and ^2 = 1-75. 

For the quiescent state we have large uncertainties on the position of the Compton peak: 
the boundaries of region A refer to values of z/c ~ 3 • 10^^ Hz (upper line) to z/c ~ 3 • 10^^ Hz. 

Since a measurement of a time lag between variations in different frequencies is lacking 
the associated constraint cannot be used. 

The longer variability timescales and the lower luminosity with respect to the two 
previous sources allow in this case a value of 6 around 10, lower than those inferred in Mrk 
421 and PKS 2155-304, as can be seen from Fig. (^, while the value found for the magnetic 
field intensity, B ~ 0.3 is similar. 




(35) 
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For the high state we must use the formulation in the KN limit: in fact the Thomson 
condition ( ]T3| ) yields a value of 5 ~ 150, highly implausible. For the observables we have 
used the following values: z/^ = 10^^, PgLsiys) = 5 ■ 10'^^ erg/s, UcLdi^c) = 1-5 ■ 10'^^ erg/s and 
the variability time and the spectral indices as in the quiescent state. The Compton peak 
frequency is not well determined since the spectra measured in the TeV range suffer from 
various uncertainties (Samuelson et al. 1998). However all the measurements give slopes 
steeper than 1 indicating that the peak is at lower energies. We therefore consider values 
between = 10^^ - 10^^ Hz. 

In the KN regime the parameter region allowed by the peak frequency ratio (region 
A) has B increasing with 6 [see Fig. (^]. The luminosity ratio (region B) has almost no 
intersection with region A for values of 6 allowed by the transparency condition. On the 
other hand the emitting region could be smaller than imposed by the variability time scale 
by a factor larger than 3. In that case region B would extend upwards defining a permitted 
region ati?>l, 5>7 — 8 

However, the condition of consistency for the value of — 10^, given in the KN 
limit by eq. ( p^ ) (region C), is definitely outside the parameter region allowed by the 
other constraints, implying that the shape of the electron spectrum cannot be determined 
balancing radiation losses and escape, unless again the size of the region is much smaller 
than estimated from the observed variability. Since the flaring state lasts much longer than 
the cooling time, efficient continuous reacceleration processes are needed. 

A comparison of the parameter space for the low and high states of Mkn 501 indicates 
that not only 7^ but also B increased in the flaring state. 

7. Summary and Conclusions 

We have extended the analytic treatment of the spectral constraints on the 
SSC homogeneous model for blazars initiated by GMD including KN corrections and 
adding the discussion of other constraints now available. The constraints can be 
summarized in the magnetic field-Doppler factor plane. In principle the knowledge of 
J^s,peak,^c,peak,Ls,peak,Lc,peak plus t^ar allows the determination of B, 6, -fb, R and K. The 
transparency condition, the time lag condition, the cooling-escape balance are further 
restrictions that overconstrain the model. 

In practice large uncertainties on the observational quantities still exist, but at least 
for the best observed TeV blazar, Mrk 421, we can significantly constrain the parameter 
space. We find that the radiative interpretation of the lags is consistent with the spectral 
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constraints. Moreover, the break in the distribution is consistent with the equihbrium 

between the energy losses and electron escape timescale. The derived value of the the 
Doppler factor is moderately large {S ~ 25), in agreement with the values found by specific 
spectral fits (see e.g. GhiseUini & Maraschi 1996; Mastichiadis & Kirk 1997). 

We have apphed the analysis also to PKS 2155-304, a well observed object in the low 
energy bands, detected in two occasions with EGRET in 7-rays. Assuming consistency 
between all the constraints (as found in the case of Mrk 421) we estimate that the peak 
of the Compton component should fall at about 10^^ Hz: for higher values we do not find 
consistency between the spectral constraints and the time lags as measured in the 1994 
observational campaign (Urry et al. 1997). We estimate the TeV flux, which could be 
measured by southern Cherenkov telescopes. 

The last object analized is Mrk 501, a TeV emitter that has shown a dramatic flare in 
April 1997. For the quiescent state we can apply the classical formulation, while for the 
flaring state we must use the formulae in the KN limit. The absence of information about 
time lags prevents to completely constrain the model: however we can estimate that the 
value of the Doppler factor is about 10, smaller than in the other two sources. The high 
state is consistent with a costant Doppler factor, but impUies an increase in ^yj, and in the 
magnetic fleld. The required value of 76 is inconsistent with cooling-escape equilibrium, 
indicating the need of an efficient continuous reacceleration mechanism. 

In conclusion the analytic approach offers a solid basis for a general discussion of the 
model and its applicability to specific sources. It allows to derive parameter ranges taking 
into account uncertainties in the "observable" quantities and a better comprehension of the 
interrelations between various parameters. We stress however that Table 2 has been derived 
applying a number of restrictive choices and hypothesis for illustrative purposes and should 
not be taken at face value. 

Individual fits remain essential to better understand the detailed shape of the particle 
spectra and their origin. 

Note: after the paper was accepted we learned about the detection of PKS 2155-304 
by the Durham Mark 6 Cerenkov Telescope (Chadwick et al. 1998) with a fiux of 
F(> 0.3re\^) = 4.2 • 10~^^ ph cm~^ s~^, consistent with our estimate. 
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A. The Klein-Nishina limit 

In this Appendix we give the expression of the Compton emissivity for the electron 
distribution assumed in this work and we prove eq. (|np for the frequency of the Compton 
peak in the KN limit. All the quantities are referred to the blob reference frame. 
As previously assumed, we consider a population of relativistic electrons with a double 
power-law distribution with indices rii and n2 and break Lorentz factor 7^: 

iV(7) = P^^"" 'l^^^' (Al) 

^ \ i^27-"^ if 7fe < 7 < 7ma. ^ ' 

where 'jmax > lb- 

The coefficients Ki and K2 are related by the condition ^^17^"^ = K2'~f^"''^ which gives 

if2 = i^i7r""^- 

The seed photons are those emitted by the synchrotron mechanism and in the blob reference 
frame the radiation energy density can be written as (hereafter we use the parameter 
x' = hu'/mc^) 

eo -a for x^^^ >Xq>x^ 



where a, = (rij — l)/2. 

Using a ^-approximation (see e.g. Coppi & Blandford 1990) the Compton emissivity in the 
blob reference frame can be written as: 



Stt/i Jxi 



3x'V^^ eix'r 



Ax'n J x: 



^dx'o (A3) 



where the limits of the integral are analyzed in the following. 

For iV(7) and e(xo) given above we can split the integral in four similar integrals: 

j{x')=Ii,l + Il,2 + l2,l + l2,2 (A4) 

The first term is: 

h,, = C{a^)K,eoK] ' In ^ (A5) 

V X, I X2 
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with C{ai) = (4/3)-°i((7Tmc^)/(87r/i). The limits are: 



Xi = max 



X. 



3x' 



; X2 = mm 



^ 3^ 3 
T' 4^ 



X, 



(A6) 



The term 3 /Ax' in the upper hmit express the KN suppression of the high energy emission. 
It originate from the step approximation of the KN cross-section, that can be rewritten as 
the condition that the scattering are in the Thomson regime up to: 



This condition, with the use of 



gives: 



7^0 < 4 



x' = {4/3hX 
. 3 1 



Xr, < 



Ax' 



that is the condition used in the integrals. 
The second integral is: 



h,2 = C{ai)KieQX' 



l—ai IOL2 



OLi — 0L\ 



ai— 02 ™Q!1— ^2 



where: 



Xi — max 



X 



, 3x' 



SI A 2 



; X2 — mm 



3x^ _3_ 



The third term can be written as: 



(A7) 

(A8) 
(A9) 



(AlO) 



(All) 



72,1 = C(a2)Xi7?°^~"^^eo^'-"^<"^^ 



0^2 ~ Oil 

where, in this case C{a2) — {A/3)^°'^{aTmc^)/{87rh) and: 



a2— ai ™a2— ai 



(A12) 



xi — max 



X, 



3x' 

miw 7~2 

I max 



; X2 — mm 



X 



, 3x' 3 



(A13) 



The last term is: 
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(/ \ -"2 
x'J Xi 



with: 



Xi — max 



3a;' 



X 



" 472 



; ^2 = mm 



3a;' 3 
4^' 4^ 



(A14) 



(A15) 



The Compton emissivity around the peak in the KN hmit is dominated by the 
contribution of and /2,i (see Fig. 2 d). One can easily see that the emissivity can be 
written as: 



X 



I \ — 



3{x') - C(«i)i^ieo - 
\ X „ 



In 



■ 2 1 



0^x0.2— Oil 

4 



q;2 — 0.1 



(A16) 



/^From this expression we can find the frequency of the Compton peak, that is the maximum 
of the function j{x')x'. Taking the derivative we can write the maximum condition as: 



In 



■ 2 1 



+ 



OI2 — Oil 1 — oil 
which gives the value of the frequency of the Compton peak: 



= 



(A17) 



x^ = 76 exp 



+ 



ai — 1 2(0:2 ~ Oil) 



(A18) 



B. Constants used in the text 



The constant ^(cti, 0:2) has the expression: 



C(q;i, 0:2) = 2-10-^ 



h \ i-'"! 



/(«1,"2) 



-|l/(l-ai) 



2-3ai 



mcy [5.7-1039 

In Table 3 we give the value of ^(cti, 0(2) for some common values of the spectral indices 



(Bl) 
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Figure Captions 



Fig. 1. — SED of Mrk 421 in two epochs (triangles and squares, data from Macomb 
et al. 1995). Spectra calculated with the homogeneous SSC model with the electron 
energy distribution N{^) = K^~^^{1 + 7/7b)"^^~"'^ are superposed to the data: both peaks 
are reasonably well described by this double power-law approximation. The models are 
calculated with the following parameters: low state: R = 2.7 ■ lO^^cm, B = 0.15G, 
K = 1.7 ■ 10^ 6 = 25, 7fe = 5.6 • 10^, ni = 2.2, = 4.5; flaring state: as in the low 
state but 7;, = 1.4 ■ 10^ iiT = 2.4 • 10^ Note that in the radio band (z/ < 10^^ Hz) the 
emission is selfabsorbed. 

Fig. 2. — Inverse Compton SSC spectra calculated using the formulae described in the 
Appendix (the frequency is expressed as x' = hv' /mc^). The solid line is the total spectrum, 
while the other curves show the contributions of electrons and photons with different 
energies (see the following). The spectra are obtained for fixed synchrotron peak frequency 
{v'g = lO^^ifz) and spectral indices (cti = 0.5, a2 — 1.5), but for different Lorentz factors 
of the electrons at the break. The dotted line shows the contribution from the electrons 
and the photons before the break (/i,i), the short dashed fine shows the contribution from 
the Ji 2 branch (low energy electrons and high energy photons), the long dashed line is the 
contribution from l2^i (high energy electrons and low energy photons) and finally the dot- 
dashed curve is the /2,2 spectrum (the contribution from the high energy photons and high 
energy electrons). In panel a) 75 = 10^ and the Compton peak is produced in the Thomson 
limit: all the contributions are present and the peak frequency is given by = (4/3)7^2/5. 
For higher values of 7^ [panel b)] the contribution from /2,2 begins to be suppressed; when 
the peak is in the KN regime [panel c) and d)] the latter contribution vanishes, while also 
the contribution from 7i^2 is strongly suppressed. In the strong KN limit, only the and 
/2,i branches are important. Note that the spectrum after the peak in the KN limit is a 
power-law with the same index that the /2,i component, ckkn — 2q;2 — cti 

Fig. 3. — Constraints on the parameter space for Mrk 421. Region A represents the peak 
frequency constraint. The region takes into account an uncertainty of a factor 3 in the 
position of the Compton peak, the upper line is calculated for Uc — 2- 10^^ Hz, the lower fine 
for Uc — 1.8 • 10^^ Hz and the intermediate line for Uc — 6 ■ 10^^ Hz. Region B is obtained 
from the peak luminosity constraint and takes into account an uncertainty of a factor 3 in 
the dimension of the region (the upper line is calculated for t'^^^ = l/St^ar)- The region C 
is obtained from the equilibrium of cooling and escape for electrons with Lorentz factor 7{,; 
the region is plotted taking into account an uncertainty of a factor 3 in the value of tyar and 
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in the value of /3esc (1-1/3)- The thick hne I is obtained from the abserved time lag between 
variations in different X-ray band. The vertical line shows the lower limit on S given by the 
transparency condition. The plot shows that consistency of all the constraints is achieved 
for (5 ~ 25 and 5 ~ 0.2 

Fig. 4. — Parameter space for PKS 2155-304. Lines and regions are labelled as in Fig. 3. 
Region A takes into account an uncertainty of a factor 3 in the Compton peak frequency 
(2.5- 102^.5- 1025). 

Fig. 5. — Parameter space for Mrk 501 (low state). For this object we do not have 
information about delays in the variability and therefore in the diagram we cannot plot 
the line 

Fig. 6. — Parameter space for Mrk 501 (high state). In this case we have used the formulae 
for the KN limit and region A has a positive slope. Note that the break energy-balance 
condition between energy losses and escape (region C) is not consistent with the other 
constraints. 
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Table 1. Input observational quantities 



Mrk 421 PKS 2155-304 Mrk 501 (l.s.) Mrk 501 (h.s.) 



Us (Hz) 3 • 10^6 IQi^ 10^^ 10^9 

(Hz) 2 ■ 10^5 - 2 • 10^6 2.5 ■ 10^4 - 2.5 • lO^^ 3 ■ 10^^ - 3 • lO^^ 10^5 - 10^6 

z/,L,(z/,) 4.7-10^^ 1.5-10^6 10^4 5-10^5 

i/cLe(z/c) 6.3 • 10^^ 1.5 ■ 10^6 5 • 10^^ 1.5 ■ 10^^ 

t^ar (h) 1 5.5 10 10 

r^is) 4200 4700 



^lag between photons of 5 keV and 0.75 keV 
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Table 2. Physical parameters estimated from the present analysis 



Mrk 421 PKS 2155-304 Mrk 501 (l.s.) Mrk 501 (h.s.) 



5 22 - 35 24 -38 8 -20 7 

B (G) 0.3 -0.2 0.3- 0.2 0.5 - 0.1 1 

76« 4 - 3 2 - 1.8 3 -4 60 

3 - 1 15 -7 9- 7 2.5 



^in units of 10^. 
''in units of 10^^ cm. 
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Table 3. Values of ^(cki, a^) for some values of the 
spectral indices 



0.5 1.25 9.13-10-^^ 

0.5 1.75 2.17 • 10-^6 

0.75 1.25 1.54 ■ 10-^^^ 

0.75 1.75 1.1 ■ 10-165 
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